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The gonihedric spin model was first introduced as the action for a dis- 
cretized tensionless string in a discretized embeding space. Afterwards was 
found that there are interesting features on the dynamical behavior of this 
model in 3 dimensions (as it was first formulated) that make us think on 
glassy spin model without inherent disorder. Extensive simulations have 
been carried out in the 3 dimensional model. In the following I will re- 
port on a work composed of two different but related parts. The first part 
is a numerical study through Monte Carlo simulations of the dynamical 
properties of the 2 dimensional version of the model (i.e. the loop model), 
which is much simpler due to the fact that it has trivial thermodynamical 
properties. The second part consists on an analytical approach of this 2 di- 
mensional loop model coupled to gravity. We solve partially the associated 
two-matrix model via a reduction to an equivalent one matrix model and 
saddle point methods with the last one-matrix model. 

1. An statistical model of loops in 2D 

The gonihedric spin model was first introduced by Savvidy in relation 
to a discretized model for a tensionless string theory [T], but very soon the 
spin model gained interest by itself. Also its extension to a self-interacting 
surfaces (k / 0) showed a very rich family of models with different kind 
of critical points and interesting dynamical properties [I][2]f3]. Extensive 
numerical and theoretical work appeared [4||5j[6j[7f and some interest about 
the glassiness of the 3-dimensional gonihedric model arised [8J|9J[10J|11J[12J. 
This is precisely the aspect of the model that has motivated us to study the 
2-dimensional version of the model. Its trivial thermodynamics motivates 
us to investigate whether this model also has glassy behavior or not. This 
would provide a toy model for glassy phenomena without inherent disorder 
on the couplings. 

* Presented at the workshop Random Geometry Krakow 2003 
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So the spin model we are going to investigate is related to a model of 
loops in 2 dimensions. We are going to compare this model and the results 
we obtained from analogous simulations in a 3 dimensional version of this 
spin model But before we go into the analysis of the model let's define 
it in terms of spin variables. 

Consider the following Hamiltonian 

Klonih = ~ K Yj a * a S + \ Yl ai(T i + ~Y~ Yl a i a 3°K a l (1-1) 

<i,j> <ij> [i,3,k,l] 

in a two dimensional lattice 1 , where < i,j > means sum over nearest neigh- 
bors, « ij » means sum over next to nearest neighbors, and [i,j, k, I] 
means over spins forming plaquettes on the lattice. 

This Hamiltonian has some odd characteristics. The most important 
of all is that the space of symmetric vacua is extremely large, in fact it is 
exponentially large with the dimension of the lattice L. In particular, the 
simultaneous flip of all the spins that belong to any set of non-crossing lines 
leaves the energy of the ferromagnetic ground state unchanged 2 This 
symmetry is even larger in the n = case where the lines can cross each 
other. This provides a very special landscape for the energy function of our 
model that in its 3 dimensional version makes the system exhibit a very 
clear glassy behavior associated to a thermodynamical phase transition [TO] . 
This is precisely the aim of this work: to determine whether or not the same 
kind of behavior can be found in 2 dimensions given that the 2D model has 
no thermodynamical phase transition. 

Let us first see the relation between this model and the loop model 
we announced. If we look at the energy of a given configuration we can 
see that due to the precise fine tuning of the couplings all the energy is 
concentrated at the bending points of the loop (surface in this 3D case) that 
is the boundary between the two different phases of the system (i.e. between 
plus and minus spins) , and that wherever there is a crossing of this surface 
with itself (or with another loop) there is another concentration of extra 
energy. So at the end of the day we can write the energy of the spin model 

1 The 3 dimensional counterpart has slightly modified couplings, the explicit form is 

ri gonih = -2K (TiCTj + - <Ti<Tj-\ — 

<i,3> <Ki,j> [i,j,k,l] 

The reason is that the number of neighbors changes from 3D to 2D, thus the couplings 
has to change too as we will argue. 

2 In the 3 dimensional case this symmetry is generated by the flip of planes rather than 
lines. 
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(or the loop model) as follows 

£ , = n 2 + 4Kn 4 (1.2) 

where ri2 is the number of bending points and the number of self-crossing 
points of the loop that separates the plus and minus spins regions. 

This is exactly the same that happens in the 3 dimensional version of the 
model. In fact the couplings are precisely chosen to exhibit these features. 
But all this has very different consequences for 2 dimensions or 3 dimensions. 
The main difference between this 2D loop model and the corresponding 3D 
surface model is that the action of the surface in 3D is proportional to the 
linear extent of the surface (see figEO and the roughness of it, while the 
action for the loop is not depending on how big it is but on how many times 
it bends (see fig^t). Thus in 2 dimensions the energy of the loop do not 




Fig. 1. Examples of configurations of the system in 3D (a) and in 2D (b). The 
energy is concentrated on the bending and crossing points (2D) and lines (3D). 

depend on its size but on the shape of it (although the energy barriers do 
eventually depend on the size of the loop). 

1.1. Thermodynamical behavior 

Let's now take a look at the thermodynamical properties of this spin 
model taking as a reference its 3 dimensional counterpart |K|[12j. Let's 
begin with the special case of At = that is exactly solvable in infinite 
volume and reducible to an easy-computable sum for finite volume. The 
exact solution for the model with k = shows us that there is no phase 
transition at finite temperature. If we take a look at figEl we will see the 
infinite volume energy function and susceptibility compared to the numerical 
results of simulations and to the exact finite volume calculation. All the 
discordances between simulations and the infinite volume calculations are 
due to finite volume effects as we can see comparing the simulations with the 
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Fig. 2. (a) Energy function and (b) specific heat of the system for k — 0. The exact 
function at infinite volume, at finite volume, and the Monte-Carlo simulations are 
plotted 



exact finite volume calculations. For the other cases with there is no 

infinite volume exact solution nor easy-computable finite- volume expression 
but the simulations performed do not show great differences with the k = 
case (see figEJ- The only remarkable difference is the appearance of a second 




Fig. 3. (a) Energy function and (b) specific heat of the system for k = 0. The exact 
function at infinite volume and the Monte-Carlo simulations are plotted 

structure for sufficiently large k. This second structure can be interpreted 
as the appearance of a new energy level for the plaquette variables. This has 
been studied to see whether it evolves into a peak at large volumes, but no 
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volume dependence of this structure has been found, so there is no evidence 
of phase transition. 

On the other hand the same model in three dimensions exhibits a quite 
complex phase space. For k = there is a critical temperature T c where 
the system changes from ordered to disordered phase through a first order 
phase transition, and a second temperature T g that is between two differ- 
ent dynamical phases: a glassy phase and a supercooled phase (see fig@J 
E] Increasing k we find from certain value on that this second tem- 




1.5 1.6 1.7 1.8 1.9 2 2.1 

Temperature 



Fig. 4. Energy versus temperature. The lower branch is produced heating the 
ferromagnetic ground state. The higher branch corresponds to sudden quenches at 
each temperatures from disordered configurations. 

perature T g either is very close to the thermodynamical temperature T c or 
they coincide, and that this thermodynamical phase transition changes from 
first to second order. 

Thus we have seen that this models from 3 dimensions to 2 dimensions 
changes a lot its thermodynamical behavior. From having first/second order 
phase transitions (depending on the value of k) on 3D to trivial thermody- 
namics without any phase transition on 2D. But we also want to know 
whether there is a great difference or not in they dynamical properties. In 
particular we want to know if the slow dynamics and the glassy behavior 
remains on the 2 dimensional model. 

1.2. Dynamical behavior 

We will move now to the dynamical properties of the model. We shall 
report here only the K = case 3 . 

3 Due that we have not established yet without ambiguities whether this case posses 
or not glassy behavior we are not going to consider the case. 
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To study the dynamics of this system we will consider a two-time corre- 
lator of local observables ^Hj 

C(t,t w ) = ^2ei(t w )ei(t w + t) (1.3) 

i 

where the sum runs over all the sites in the lattice, and the variable ej(i) is 
the energy 4 of the site i at the time t. This object, in equilibrium, should 
be independent of t w , but as we can see in fig El below some temperature 




(a) (b) 



Fig. 5. Autocorrelation functions for different temperatures and different waiting 
times: (a) curves are independent at t w (b) some dependence in t w appears 

this function happens to depend on the waiting time t w , such observation is 
identical to the one made for the 3 dimensional case where an even larger 
dependence of this autocorrelation function on the waiting time appeared 
below T g (see figEI- This is not per se a clear evidence for glassy behavior 
of the model yet, so we will continue with the program we followed with 
the 3 dimensional model. Thus let's fit the curves that look i^-independent 
(i.e. the curves that should be above the hypothetical glassy transition 
temperatureTg). The fitting function will be of the form 

Ae&" (1.4) 

which is an stretched exponential. As you can see in figEK the agreement 
between the fit and the simulations is rather good 5 . Now we can plot the 

4 We could have used other kind of observables like the spin variables or the energy 
per plaquette, but they have the same behavior for our purposes. 

5 Lines are fits, points are simulation measurements. We looked for the consistency of 
the fits by checking that the true value of the parameter A, i.e. A = 1 (which we know 
by construction of the correlation function) were within the interval of confidence of 
the fitted value 




Fig. 6. curves of autocorrelation of the spin variables at both sides of T g on the 3 
dimensional gonihedric model 



fitted values of r against temperature, and we see that it increases as we 
lower the temperature as if it liked to diverge at some point. If we fit this 
points using a function of the form 

{t-°%y (L5) 

as we did in 3 dimensions, we'll find a good parameterization of the di- 
vergence (see figd although there is still something that is not completely 
compatible with this fit. The problem with this fit is that the fitted "critical" 

10 I n 1 1 1 1 j 1 1 

\ Measured values i — I — i 
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Fig. 7. Autocorrelation time versus Temperature at temperatures above 0.9 
temperature T* is not consistent with the point where the autocorrelation 
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function began to depend on the waiting time, as it should. This is why we 
should analyze more carefully this autocorrelation function. Looking more 
carefully at this autocorrelation function we can see that the dependence in 
the waiting time t w disappears as we increase it. This means that the ther- 
malization of these two time functions is extremely slow, but that could be 
non- glassy-like. In fig|E]we can see this convergence of the autocorrelation 
functions as we increase the waiting time. 



C(t,100) 

C(t,1000 

C(t,10000) 




Fig. 8. Evolution of the autocorrelation function with the waiting time t w . 

To explore better the dynamics behavior of this model we can perform 
other type of experiments in our system. For example we can prepare our 
system in a specific configuration, for example a lattice with two different 
coexisting vacua, one inside the other (one possibility could be layer-like 
vacua inside ferromagnetic one) and look at the decay process of the system 
to the equilibrium. The problem with this tests is that as there is no ordered 
phase in this model, we cannot prepare the system in an initial configuration 
composed by two different coexisting vacuas and pretend that they are more 
disordered than the equilibrium-like configurations at the temperatures we 
are examining, as it is in the 3 dimensional case (figEl shows the 3D decay 
of a perturbed vacua to the unperturbed one. It is easy to see the different 
behaviors in terms of temperature), in spite of this we performed those 
simulations and found that the decay behaved in the same way for all the 
range of temperatures ffiglTUh and figllOb are two examples of this decays 
at both sides of the hypothetical T g ). The magnitude we used to study this 
decay in two dimensions is the energy difference with the equilibrium. If 
we look at the value of the fitted exponent c of the function l|1.4jl that we 
have also used in this experiments, we can see that they are really close to 
one, and this suggests that the behavior may not be glassy but 'usual' albeit 
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Fig. 9. In the 3D gonihedric model: (a) Evolution of some order parameter start- 
ing from a perturbed vacua as an initial configuration. We can see two kind of 
behaviors. In (b) we see that the low temperature behavior is logarithmic. Lines 
are evolutions for different initial volumes of the perturbations at a temperature 
deep inside low temperature region. 




Fig. 10. The decay of a prepared initial configuration to the equilibrium 



rather slow exponential decay. 

In the analysis of the glassiness of the system we can introduce a new 
observable; the Q parameter This parameter is defined in the following 
way: after evolving a single system t w Monte Carlo steps, we make two copies 
of the system and let them evolve independently, then the Q of a local 
observable is the overlap of this observable between the two independent 
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copies of the system. In our case we will use the local spin variables 

Q = ^2^{t w + t)a b i (t w + t) (1.6) 

i 

Then, using this Q parameter and the time overlap of the same local mag- 
nitude, 

C S pi n (t : t w ) = ^ &i(t w )(Ji{t w + t) (1-7) 

% 

we can perform different kind of analysis. One of those is the following. It's 



C It, 100) 
Qjt/2, 100 

c(t, iooo 

Q(t/2, 1000 




C(t, 100) 
Qlt/2, 100 

c<t, iooo 

Q(t/2,1000 
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(b) 





(0 



(d) 



Fig. 11. Plots of the functions C(t,t w ) and Q(t/2,t w ) for different temperatures 
and different t w . In (a) and (b) we see how the scaling 11.80 is satisfied for any t w . 
In (c) and (d) we see that the relation of scaling is not satisfied at low temperatures. 

known that the system must satisfy the relation 

C S pi n (t w ,t) = Q(t w ,t/2) (1-8) 
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if the we are in an ordinary non-glassy phase. As we can see in fig lTTh , and 
fig lllb there are temperatures 6 for which the behavior of the Q parameter 
with respect to the C sp i n two times auto-correlator is what we expected, 
while there are lower temperatures, like in fig lllb and fig lllH where the 
relation (| 1 . H|) is not satisfied at large enough times, for the value of t w we 
simulate, although the discrepancy region is moving to larger times as we 
increase t w as if it would like to follow (jl,8|) for large t w (see figll2jh 



Fig. 12. The convergence to the scaling behavior is clear in this plot where t w has 
been increased by one order of magnitude. 



At this point we have concluded the analysis of all our simulations with- 
out any clear reason to believe that this 2 dimensional version of Savvidy's 
gonihedric model has really glassy behavior; on the contrary it seems to 
exhibit only very slow dynamics not related to real glassiness of the model. 
This has to be further investigated to clear out what type of dynamics is 
this model developing. Also has to be investigated although we think 

it will follow the same kind of behavior that the k = case. 



There exist a way to extent this model to one coupled to gravity. To 
couple it to gravity we will put our spin model in a random lattice built 
from quadrangular pieces. In this way we will keep the behavior of the 
loops and add the gravity degrees of freedom. In order to make the matrix 

6 Remember that T = 0.8 is already a temperature where, at waiting times we are 
considering, the energy auto-correlators are not t m independent 



0.01 




t(MC steps) 



1.3. Conclusions 



2. The model of loops coupled to gravity 
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model solvable (or approximately solvable) we will make the loops highly 
self-interacting, that means that the loops will never cross themselves. This 




Fig. 13. Example of a random lattice with a gonihedric spin model on it 

corresponds to the k — > oo limit in the model we presented above. In figlTTfl 
we can see an example of this kind of quadr angulations. From this picture we 
can extract the weights of each interaction term in the matrix model that 
will represent the partition function of our system. Let's see what those 
terms mean. 

First of all there will be the bulk term (in other words, plaquette that is 
not crossed by any loop), then we have to consider a term were a plaquette 
is crossed by one loop without bending through it and finally the term were 
the loop crossing the plaquette bends in one direction or the other. These 
three building blocks of the random lattices are presented graphically in 
fig ll4l with the corresponding term of the matrix model that will generate 
them 7 . As we can see we are considering the most general case where all the 



gABAB 



X'A~ B~ 



g'A 4 



Fig. 14. correspondence between the loop pieces and the matrix interaction that 
are going to generate them. 

couplings are different but in our specific case we will impose the condition 
g = g due to the fact that a straight piece of loop do not contribute with 
any amount of energy to the action, so the coupling has to be equal to the 
bulk coupling. In addition to those terms we have to put the kinetic term 

7 To simplify the visual identification with the loop model we have not drawn the lines 
corresponding to the "bulk" propagator, i.e. to the A matrix propagator. The loop 
is generated with the B matrix propagator. 
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for the two matrices, i.e. the quadratic terms A 2 and B 2 . So finally the 
matrix model that will represent our loop model coupled to gravity will be, 



J j dALB exp( - Tr[A 2 + B 2 + j^A 4 + ^A 2 B 2 + LaBAB]) (2.1) 



Now that we have found the matrix model that will reproduce our loop 
model coupled to gravity we only need to develop its solution. To our knowl- 
edge this model has not been solved exactly. Although very similar matrix 
models have been indeed solved ^H] their solution cannot be applied to our 
matrix model. 

2.1. A partial solution to the matrix model 
Let's proceed then to the approximation to the solution. To this aim we 
are going to rescale the A matrix in the form A — > \J~^fA so that the action 
will have the following appearance. 

S = ^ Tr[yl 2 + A 4 ] + Tr[B 2 + XA 2 B 2 + gABAB] (2.2) 

where we have made some redefinitions like X'/g — > A and g'/g — > g. Now 
we are going to pay attention to the -B-dependent part. As the action S is 
quadratic in B we should be able to integrate out the B matrix and find a 
one matrix model equivalent to the one we are using now. 
The integration we are facing now is the following 

J dB exp( - Tr[£(I + XA 2 )B] - g Tx[ABAB}) (2.3) 

Here we are going to interpret the first part of the action as the free action 
(and rename I + XA 2 = M) and the second part as the interaction, so we 
can do perturbation theory and re-sum all the diagrams at the end. But 
before we calculate diagrams we need to know the free propagator of the B 
matrix, and to reach this we add some external currents and perform the 
quadratic integration 8 . So finally we find the propagator 



(BijBkt) = Z(0, AO 



1(g)! 



I <g> M + M <g> I 



(2.4) 

il;kj 



where Z(0, M) a determinant coming from the B integration. Once we have 
found the propagator we can proceed with the diagrammatic. We will only 
consider the connected diagrams. The Feynman rules for the diagrams will 
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Ax A 



B B 



Fig. 15. Feynman rules for our matrix model 

be those shown in fig ll5I To order n in the interaction term there will be also 
a factor —g/nl due to the expansion of the exponential, and a combinatorial 
factor of (n — 1)! 2 n ~ 1 coming from the reordering of the interaction terms 
( in figE3we can see the kind of diagrams that will contribute). So finally 




Fig. 16. two examples of connected diagrams that contribute to the integral, a 
non trivial contraction is shown in the second diagram, can be easily seen that is 
exactly equivalent to the trivial diagram. 

all connected diagrams add up to 



E 

n=l 



(-g) n (n- 1)! 2™" 1 



Tr 



n 

oo 

E 

n=l 



n 



2 



log[l<S)I + 



Tr (A®A[I®M + M® I] -1 )" 

2gA <g> A 
2I&I + A 2 + A 2 & I) 



In appendix^]is shown in detail how to make this calculation and find Z(J, M). 
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where we have recovered the explicit form of M in terms of A. Thus ex- 
ponentiating this last expression we recover all connected and disconnected 
diagrams; i.e. the integral 1)2.3)1 we were trying to calculate. 

J dB exp{ - Tr[B(I + \A 2 )B] - g Tr[ABAB]) = 

2gA ® A 



Z(0,M) expl — Tr 



log 1(g) I + 



21 <g> I + A(I <g> A 2 + A 2 (g> I) 



So finally we have found an expression for the integration of the B matrix. 

„ _ AT 

In this expression Z(0, M) = [det(I + XA )] 2 that can be included in the 
effective action as — y Tr[Zog(I + A^4 2 )]. This means that we have rewritten 
our two matrix model in terms of a one matrix model. So now we can 
diagonalize our remaining A matrix and integrate over the rotational degrees 
of freedom to leave our partition function in the form 



dAdB exp( - Tr[A 2 + B + g' A A + XA B + g'ABAB 



I CjJ d ai )A 2 (a)exp{-^S eff ({a})} 
J \ i=1 J 9 



where A (a) = rL<j( a i 
S e // is the following 



a,) is the usual Van der Monde determinant and 



iY 



Se//(W) = E ( a * + < + f Ml + Aa 4 2 )) + 

2gajaj 



i=l 
AT _ 



2 + A(o? + a?) 



Since we are working in terms of eigenvalues we can use standard proce- 
dures to try to solve the model. Using saddle point approximation we arrive 
to a set of coupled equations that look quite difficult to solve exactly. In the 
N — ► oo limit those equations read 



a + 4tT + g 



A a 



1 + Aa 2 



+ 



dbp(b)- 



2gb(2 + X(b 2 -a 2 )) 



2g 



idb^L 
J a-b 



(2 + X(b 2 + a 2 ) + 2gba)(2 + X(b 2 + a 2 )) 
g(tu(a + ie) + u>(a — ie)) , 



where as usual p(a) = UrriN^oo jf^i=i^( a ~ a i) an( l the resolvent uj{z) is 
defined to be equal to J db p(b)/(b — z) This self-consistent equation has to 
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be analyzed carefully to see whether there is any fixed point in the set of 
parameters that allowed us to pas to the continuum. 



In the conclusions of this second part we may comment that, since the 
fixed geometry model we have studied in the first part does not posses any 
thermodynamical singularity, even for k — > oo, we would naively expect the 
present matrix model not to exhibit any scaling limit, but this issue deserves 
further analysis. Related to this matrix model there are other matrix models 
that can be exactly solved ^E]. Although those other models have some 
critical differences, their solutions may give some hints on how to exactly 
solve the model we are interested in. In fact some of the solved models can 
be found as an special limit of ours. That could be used as a test or a guide 
to find the solution. 



We express our gratitude to the organizers of the workshop on Random 
Geometry Krakow 2003 for a most enjoyable atmosphere and hospitality. I 
would like to thank A.Dominguez and L. Tagliacozzo for interesting discus- 
sions and encouragement. The research contained in this talk is supported 
by the European Network EUROGRID and for a CIRIT grant 2001FI-00387. 



To find the propagator we begin with the free action and add an external 
source. In our case this will lead to 



Then, to reabsorb the external field J into the B field we do a linear change 
of variables B — > B + C and quadratize the action. Doing the inverse 
procedure we find 



2.2. Conclusions 



Aknowledgments 



Appendix A 

Finding the B matrix propagator 



TV [BMB + JB] 



(A.l) 



BMB - CMC ={B- C)M(B - C) - CMC 
= BMB - (CMB + BMC) = BMB - JB 



So at the end the integral is 



Z( J, M) 



dB exp{ - Tv[BMB - JB}} 




(A.2) 
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where the C matrix can be determined from the equation 

CM + MC = J 



(A.3) 



Let's solve this equation to find the explicit form of the C matrix. To this 
aim let's choose the basis where the M matrix is diagonal. So if f2 is the 
matrix that diagonalizes it 



M = VLDO) where Djj = 5ijdi 
J = nj'tf 

c = nc'rt 



(A.4) 



The primed matrices correspond to the non-primed ones after the rotation. 

So now we will continue solving the equation l|A.3|l by writing it in the 
M-diagonal form, 

(CD + DC% = C[Aj + diCi, = jL 



C'ij 



V 



di + dj 



Now introduce it in eq. ([A.2|) and invert eq. ()A.4|) to find 



-i 



N 



Z(J, M) = \ det[M]j exp<^ ^ 

d; 



J ' j -d, ^ 



di + dj 3 dj + di 



N 



det[M]\ *expi ^ 



' x (di + dj) 



2 i^ik J k Aj] P™4^ 



That is the last expression we'll write for Z(J, M). Prom here we can de- 
duce all necessary correlators, like the one we are looking for; the propagator 
(BijBki) = SZ(J, M)/5Jji5J[k \ j=o- Calculating those variations and rotat- 
ing back the D matrices we find 



Z(0,M) 



N A 
m,n=l " L n 



Z(0,M) 
Z(0,M) 



M <8>I 



[I ® M + M <g> I] 2 



+ 



I® Af 



[1(g) M + M <g>If 



il;kj 



I <g> M + M <g> I 



(A.5) 
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